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ABSTRACT 

With a view to understanding the formation of double neutron-star binaries, we in- 
vestigate the late stages of evolution of helium stars with masses of 2.8 - 6.4 Mq in 
■ binary systems with a 1.4 Mq neutron-star companion. We found that mass transfer 

\jQ ' from 2.8 - 3.3 Mq helium stars (originating from main-sequence stars with masses 

\Q , of 10 - 12 Mq which underwent case B evolution, or 9 - 10 Mq which experienced 

' case C mass transfer) as well as from 3.3 - 3.8 Mq in very close orbits (P orb <0 d 25) 

will end up in a common-envelope and spiral-in phase due to the development of a 
convective helium envelope at the end of the calculation. If the neutron star has suf- 
] ficient time to complete the spiraling-in process in the envelope of the helium star 

before the core collapses, the system will produce very tight double neutron-star bi- 
naries (P rb ~ d 01) with a very short merger timescale, i.e. of the order of 1 Myr 
or less. These systems would have important consequences for the detection rate of 
Q ■ gravitational- wave radiation and for the understanding of 7-ray burst progenitors. On 

?— 1 I the other hand, if the time left until the explosion is shorter than the orbital-decay 

c/2 , timescale, the system will undergo a supernova explosion during the common-envelope 

phase. 

Helium stars with masses 3.3 - 3.8 Mq in wider orbits (P rb ^ d 25) and those 
more massive than 3.8 Mq do not develop a convective envelope and therefore are not 
expected to go through common-envelope evolution. The remnants of these massive 
helium stars are double neutron-star pulsars with periods in the range of d l — l d . 
This suggests that this range of mass (originating from main-sequence stars more 
massive than 12 Mq which underwent case B evolution, or more massive than 10 Mq 
which experienced case C mass transfer) includes the progenitors of the galactic double 
neutron-star pulsars with close orbits (B1913+16 and B1534+12). A minimum kick 
velocity of 70 kms -1 and kms -1 (for B1913+16 and B1534+12, respectively) must 
have been imparted at the birth of the pulsar's companion. The double neutron stars 
with wider orbits (J1518+4904 and probably J1811-1736) are produced from helium 
star-neutron star binaries which avoid Roche-lobe overflow, with the helium star more 
massive than 2.5 Mq, i.e. the remnants of main-sequence stars more massive than 
10 Mq in relatively wide orbits. For these systems the minimum kick velocities are 
50 kms" 1 and 10 kms" 1 (for J1518+4904 and J1811-1736, respectively). 

Key words: stars: evolution - binaries: general - stars: neutron - pulsars: individual: 
B1913+16, B1534+12, J1518+4904, J1811-1736 - supernovae: general 
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1 INTRODUCTION 

In a work published earlier (Dewi et al. 2002, hereafter called 
Paper I), we discussed two types of mass-transfer evolu- 
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tion from a helium star to a neutron star, i.e. case BA (in 
which the onset of mass transfer occurs during helium core 
burning) and case BB (the star fills its Roche lobe after he- 
lium core burning is terminated, but before the ignition of 
carbon). We found that dynamically stable case BA mass 
transfer can take place from helium stars less massive than 
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3.0 Mq. The remnants of this type of evolution are heavy 
CO white dwarfs. Case BB mass transfer from helium stars 
less massive than 2.6 - 2.8 Mq produces white dwarfs (CO 
white dwarf for masses < 1.9 — 2.0 Mq/ helium stars or 
ONeMg white dwarf for larger masses), while more mas- 
sive helium stars leave neutron stars as their remnants 1 . 
In Paper I we suggested that the systems with relatively 
high- mass helium stars and/or wide orbits are progenitors of 
type lb supernova, and lower-mass helium stars or systems 
in close orbits produce type Ic supernovae. We were also 
able to constrain the possible progenitors of the observed 
galactic double neutron star (DNS) pulsars B1913+16 and 
B1534+12, if we assume that these DNSs were produced 
from helium star-neutron star binaries which went through 
a mass-transfer phase. 

The calculations in Paper I were done up to various 
stages of evolution, i.e. ranging from the appearance of the 
first carbon-burning convective shell to the outward penetra- 
tion of the helium burning convective shell into the helium 
envelope. Although we could not follow the evolution fur- 
ther, we are interested to investigate the possible outcomes 
that the systems might produce, e.g. whether a common en- 
velope (CE) and spiral-in will occur, whether the system will 
survive the CE phase, whether the supernova (SN) explosion 
will take place after or before the neutron star completes the 
spiraling-in process, and which type of SN the system will 
produce. We also want to find the possible progenitors of the 
wide-orbit DNS J1518+4904 and to put more constraints on 
the progenitors of B1913+16 and B1534+12. These are the 
problems we will discuss in this paper. 

As the completion to Paper I we first evolved helium 
stars in wider binary systems with a 1.4 Mq neutron-star 
companion, such that Roche-lobe overflow (RLOF) is initi- 
ated during carbon core burning or beyond (which we call 
case BC mass transfer). A study of this particular type of 
evolution has been done by Habets (1986) who evolved a 
2.5 Mq helium star with a 17 Mq main-sequence compan- 
ion in a period of 20 ( ?25. Because we are mainly interested to 
study the systems which will produce DNS, we limit our cal- 
culations to helium stars in the mass range of 2.8 - 6.4 Mq. 
Helium stars more massive than 6.4 Mq do not expand very 
much, so that their evolution in wide orbits will be similar 
to the evolution of a single helium star. 

We discuss the results of case BC evolution in Sect. |2] 
the possible remnants of case BB and case BC evolution in 
Sect. |21 and the formation of the observed galactic DNSs in 
Sect. 0] The conclusions are given in Sect. |S] Throughout 
the paper, in order to avoid confusion about the various 
stages of evolution, a subscript i is used to indicate the initial 
parameters (at the start of calculation), o the parameters at 
the end of the mass-transfer phase (at the initiation of the 
CE phase) , t those at the end of the spiral-in phase (before 
SN explosion) and / the post-SN parameters. 



2 RESULTS: CASE BC MASS TRANSFER 

2.1 A brief description on the method of 
calculation 

We used the Eggleton code (Pols et al. 1995 and references 
therein) for the evolution of helium stars, assuming a chem- 
ical composition of (Y = 0.98, Z = 0.02) and without en- 
hanced mixing (the STD model in Pols (in preparation)). 

The non-conservative orbital evolution of the system 
is assumed to be affected by gravitational-wave radiation 
(Landau & Lifshitz 1958) and by the loss of mass with an- 
gular momentum from the system (van den Heuvel 1994; 
Soberman, Phinney & van den Heuvel 1997), such that 

d 64G 3 M hbMnsMt 2pq 2 -q + (a- 1)] M Ho 
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where MHe, Mns, and Mt are the masses of the helium 
star and the neutron star, and the total mass of the system, 
respectively, q = Mhc/Mns is the mass ratio of the system, 
a is the orbital separation, G is the constant of gravity, and 
c the speed of light in vacuum, a is the fraction of mass 
lost from the helium star in the form of fast isotropic wind, 
and f3 is the fraction of mass ejected isotropically from the 
vicinity of the neutron star. The stellar wind mass loss is 
given by eq. (2) in Wellstein & Langer (1999) multiplied by 
a factor of 0.5, i.e. 
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in Mgyr -1 . We assumed that during the detached phase, 
the neutron star does not accrete matter from the stellar 
wind. During the mass transfer phase, the transferred matter 
is accreted up to the Eddington limit for helium accretion 
(3 x 10 -8 Mq yr _1 ); the rest is lost from the system with 
the specific orbital angular momentum of the neutron star. 

We devide case BC mass transfer into two mass ranges, 
i.e. (i) 2.8 - 3.2 Mq helium stars, in which a partially de- 
generate core develops at the end of the calculation, and 
mass transfer becomes dynamically unstable due to the de- 
velopment of a convective helium envelope, and (ii) 3.6 - 
6.4 Mq helium stars, in which the core is only weakly de- 
generate and a neon burning convective core develops at the 
end of the calculation. Typical examples for the two mass 
ranges are presented in Figs. 0-01 and the overall results 
are presented in Tabled Until the end of the calculation the 
neutron star does not accrete a significant amount of mat- 
ter, its mass remains 1.4 Mq and therefore is not displayed 
in the table. Detailed discussions of the evolution of helium 
stars in binary systems have already been given in Paper I, 
and therefore we will focus our discussion here on the late 
stages of the evolution. 



1 These limiting masses depend on the initial orbital period of 
the system. 



2.2 Roche-lobe overflow from 2.8 - 3.2 Mq 
helium stars 

In all systems considered, mass transfer takes place on the 
thermal timescale. In contrast to case BB evolution, the 
mass-transfer rate is lower at larger initial period (see Ta- 
ble Q. This is probably because mass transfer is initiated 
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Table 1. The binary parameters in case BC mass transfer from 2.8 - 6.4 Mq helium stars. The columns give: the initial mass and 
period; the duration, amount of mass removed from the helium star, and the maximum mass-loss rate during the Roche-lobe overflow; 
the mass and period at the end of the calculations, the final mass of helium in the envelope; and the final CO and ONeMg core masses. 
Masses are in solar units, the mass-loss rate in Mqjt -1 , time in yr, and periods in days. The number in the Note column indicates 
during which consecutive convective carbon shell-burning phase the calculation stops: a implies the appearance of the shell, b during the 
shell burning, and d the disappearance of the shell. 
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3.6 0.6 3.32 x 10 3 0.554 3.7 x 10~ 4 

1.2 8.67 x 10 2 0.036 7.6 x 10~ 5 

4.0 0.5 2.85 x 10 3 0.447 3.6 x 10" 4 

1.0 6.35 x 10 2 0.020 5.4 x 10" 5 

5.0 0.4 1.78 x 10 3 0.133 1.3 x 10~ 4 

0.6 7.92 x 10 2 0.029 7.4 x 10" 5 

6.4 0.2 1.34 x 10 4 0.052 1.1 x 10" 4 
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at a very late stage and there is insufficient time to reach 
the full thermal-timescale rate. Fig. shows the evolution 
of a 2.8 Mq helium star with a neutron star in a 10 d orbital 
period. A convective shell develops in the helium envelope 
before the third carbon-burning convective shell appears. 
The helium-burning convective shell causes a decrease of 
the mass of the CO core, Mco (defined as the central mass 
with helium abundance less than 0.1, which coincides with 
the lower limit of the helium-burning convective shell in the 
envelope). The convective shell penetrates outward both in 
mass and in radius (at the end of calculation, the upper 
boundary of the convective shell is located at 98 per cent 
of the total mass and 54 per cent of the total radius of the 
star). Because the mass of the envelope is dominated by the 
convective shell, the helium envelope behaves as a convective 
envelope. As a result of mass loss from a star with a convec- 
tive envelope, the radius does not shrink as in the case of a 
radiative envelope. This response of the star drives an enor- 
mous increase in M, as can be seen in Fig.Q Mass transfer 
appears to become dynamically unstable which would lead 
to a CE and spiral-in phase. 

The start of each convective carbon burning phase (both 
in the core and in subsequent shells) causes the core to ex- 
pand. As a result of the mirror principle, the envelope tends 
to contract and the star shrinks. In between the convective 
phases, the core contracts and the envelope expands again. 
This is reflected in the behaviour of the mass-transfer rate, 
which is sometimes interupted by a detached phase. In sys- 
tems with closer orbits, a detached phase does not always 
occur (although the behaviour of expansion and contraction 
of the star is also noticed). The helium envelope in these 
systems is thinner such that the outward penetration of the 
helium-burning convective shell takes place relatively fast. 
However, whether a detached phase occurs or not does not 



affect the final situation, since all systems in this mass range 
are expected to end up in a CE phase. 

The evolution of the central density and temperature of 
2.8 and 3.6 Mq helium stars is presented in Fig. |5] At the 
end of the evolution of a 2.8 Mq helium star, the ONeMg 
core becomes partially degenerate. As a consequence of neu- 
trino losses, T c remains relatively constant as the core be- 
comes denser at the onset of the fourth consecutive carbon 
shell burning (marked by a circle in Fig. 0, while the shell 
with maximum temperature moves outward. We will discuss 
the final stage of evolution of helium stars in this range of 
mass in Sect. 13.21 

2.3 Roche-lobe overflow from 3.6 - 6.4 Mq 
helium stars 

Fig. [3] shows the evolution of a 3.6 Mq helium star with 
a neutron star in a 0^6 orbital period. A convective shell 
also appears in the helium envelope of stars in this range 
of mass. However, we do not find that the shell penetrates 
outward nor that M increases enormously like in the lower- 
mass helium stars (see the fourth and sixth panels of Fig.yJ. 
Although we could not follow the evolution further than 
neon ignition, we expect that a CE phase probably does not 
occur in this range of mass. 

The core becomes only weakly degenerate. Before the 
third carbon-burning convective shell appears, the shell with 
maximum temperature moves to the centre. When the cen- 
tral temperature is ~ 1.3 x 10 9 K, neon burning occurs in a 
convective core. Before the convective ignition, the central 
neon abundance has decreased by ~ 3 per cent by radiative 
burning (see Fig. ■ Convective neon core burning is found 
only in the 3.6 Mq model with orbital period of 0?6 and the 
4.0 Mq model with orbital period of 0?5. In other models, 



4 Dewi & Pols 



10.5 

10.4 

~^ 10.3 
n 

b 10.2 

Q_ 

10.1 
0.8 
0.6 
x" 0.4 
0.2 



1 days 






log(t f - t) yr 

Figure 1. The evolution of a 2.8 Mq helium star with a 1.4 Mq 
neutron star companion with an orbital period of 10 d starting 
from carbon core burning. The orbital evolution is presented in 
the first panel. The second panel shows the central abundances: 
dashed-dotted-, dashed-, solid, and dotted-lines are C, O, Ne, and 
Mg abundances, respectively. The third panel gives the stellar lu- 
minosity (dash-dotted-), and the contributions of helium burning 
(solid-), carbon burning (dotted-), and neutrino losses (dashed- 
line). The fourth panel presents the mass-loss rate; while the fifth 
panel the stellar (solid-) and Roche (dotted-line) radii in solar 
units. The sixth panel shows the evolution of the stellar interior. 
The upper and lower dotted lines are the CO and ONeMg core 
masses, respectively, and the dark- and light-shaded areas mark 
the convective and semiconvective burning regions. 



the decrease in neon abundance is also noticed, but at the 
end of our calculations the central temperature is not high 
enough yet for convective neon ignition. The onset of con- 
vective neon ignition is represented by point F in Fig. [5] The 
evolutionary track at this point is similar to point A, i.e. the 
appearance of the convective core and the rapid increase of 
the neon burning rate cause a decrease in density while the 
temperature increases. The mass of the ONeMg core at the 
ignition of neon is 1.6 Mq (in a 3.6 Mq) - 1.7 Mq (in a 
4.0 Mq helium star). Here the ONeMg core is defined as 
the central mass that contains less than 10 per cent carbon, 
which coincides with the upper limit of the carbon burning 
convective shell. 




5 6 
log P c (g/cm 3 ) 

Figure 2. The evolutionary tracks in the central density-central 
temperature plane of a 2.8 Mq and a 3.6 Mq helium star which 
go through case BC mass transfer. Point A corresponds to carbon 
core burning, B (and C): the exhaustion (and depletion) of carbon 
in the core, D: the start of convective carbon shell burning, E: the 
end of convective carbon shell burning, and F: neon ignition. The 
arrows indicate the subsequent convective carbon shell burnings. 
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Figure 3. Same as Fig.fflfor a 3.6 Mq helium star with an orbital 
period of 0'?6. 
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Figure 4. The central temperature and central neon abundance 
in 3.6 Mq and 4.0 helium stars around the moment of neon 
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ignition. Solid circles mark the ignition of neon in the convective 
core. 



similar way. Assuming that a penetration of the convective 
shell in the helium envelope accompanied by an enormous 
increase in M is the indication of the start of dynamically 
unstable mass transfer, we suggest that a CE phase also 
occurs in 3.2 - 3.8 Mq helium stars in very close orbits. 
Wider-orbit systems in this range of mass as well as he- 
lium stars more massive than 3.8 Mq probably will not go 
through a CE phase. 

A similar study of the evolution of helium star-neutron 
star binaries has recently been carried out by Ivanova et 
al. (2002). Although the latter work is similar to ours, they 
come to different conclusions. They do not find that lower- 
mass helium stars develop a convective envelope at the end 
of their evolution, even though their calculations were pur- 
portedly done up to a more advanced stage of evolution (i.e. 
up to oxygen burning) than our calculations, which extend 
up to radiative neon ignition. On the other hand, they find 
that helium stars with masses between 3.3 and 5 Mq in 
orbital periods less than 0"?3 (which they called case CEB) 
undergo stable mass transfer, but the mass-transfer rate ex- 
ceeds a critical value which is a function of the orbital pe- 
riod, such that a CE may form which probably leads to a 
merger. However, it is difficult to compare our models with 
their calculations as they do not provide the detailed interior 
structure of their models. 



3 THE POSSIBLE REMNANTS OF CASE BB 
AND CASE BC EVOLUTION 

3.1 The late stages of case BC mass transfer and 
comparison with case BB evolution 

We have discussed in Sect. l2~2l that a convective envelope 
develops in helium stars of 2.8 - 3.2 Mq which go through 
case BC evolution, and therefore they will undergo a CE 
phase at the end of their evolution. Such a convective shell 
in the helium envelope was also found in some of our case 
BB calculations (Paper I). We could not always follow the 
evolution far enough to see if outward penetration of this 
shell does take place. However, we found that the penetra- 
tion and enormous increase in M occur in helium stars of 2.8 
- 3.2 Mq with various initial periods of case BB evolution. 
Therefore, we suggest that a CE phase is a typical charac- 
teristic of the late stage of evolution of helium star-neutron 
star binaries in this range of mass. 

In helium stars more massive than 3.2 Mq, the outward 
penetration of the helium-burning convective shell is only 
found in 3.6 - 3.8 Mq stars with very close orbits (Pi < 0^25, 
i.e. case BB mass transfer) 2 . We could not draw a conclu- 
sion about the 3.4 Mq models because all calculations for 
this mass stopped before the appearance of the convective 
shell in the helium envelope, but they probably behave in a 

2 In Paper I, it was mentioned that the penetration of the helium- 
burning convective shell into the envelope takes place in helium 
stars of various masses (cf . table 4) . In some cases the appearance 
of the convective burning shell was noticed, and it was assumed 
that once the shell appears, it will penetrate outward. Later while 
working on this paper, we found that the appearance of a convec- 
tive burning shell in the helium envelope is not always followed 
by an outward penetration. 



3.2 Final stages of evolution of 2.8 - 3.2 Mq 
helium stars 

We have discussed in Sect. 12.31 that before 3.6 Mq and 
4.0 Mq helium stars ignite neon convectively in the cen- 
tre, the central neon abundance decreases by ~ 3 per cent. 
The decrease of the neon abundance in the shell with maxi- 
mum temperature is also found in 2.8 - 3.2 Mq helium stars. 
Fig. EJcompares the conditions in the centre and in the shell 
with maximum temperature for 2.8 Mq and 3.2 Mq stars. 
The maximum temperature increases with time. In 2.8 - 
2.9 Mq helium stars, the shell with maximum temperature 
moves outward (see panel (a) in Fig. El and at the end of the 
calculation the temperature is still lower than 1.3 x 10 9 K, i.e. 
insufficient to enable convective neon burning. In a 3.2 Mq 
helium star, the shell moves inward and at the end of the 
calculation the maximum temperature reaches 1.4 x 10 9 K, 
but an off-centre neon burning convective shell has not yet 
developed. Panel (c) of Fig. shows that in both cases, the 
mass of the ONeMg core at the end of the calculation is large 
enough to ignite neon, i.e. higher than 1.37 Mq (Nomoto 
1984). Hence, we argue that at the end of our 2.8 - 3.2 Mq 
calculations the stars are close to off-centre neon ignition. 
For comparison, Nomoto & Hashimoto (1988) found that a 

3.3 Mq helium core ignites neon in the centre while 2.8 - 
3.2 Mq helium cores ignite neon off-centre. They conclude 
that 3.3 Mq is the critical mass between off-centre and cen- 
tral ignition of neon. 

After off-centre neon ignition, the future of a partially- 
degenerate ONeMg core depends on whether the burning 
shell is able to reach the centre (Nomoto & Hashimoto 1988). 
If it does reach the centre, then subsequent nuclear burning 
stages of O and Si will form an iron core, which will collapse 
due to photodisintegration. Otherwise, part of the ONeMg 
core is left unburned and becomes highly degenerate, and 
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Figure 5. The location of the shell with maximum temperature (a), the density of this shell (6); the corresponding mass of the ONeMg 
core (c), and the central temperature (d) for 2.8 Mq (solid line) and 3.2 Mq (dotted line) helium stars, presented as a function of 
maximum temperature. 



the core will eventually collapse due to electron captures. 
The density at the burning front also plays an important 
role in determining the future of the core. If the density is 
higher than 10 8 gcm -3 , neon shell burning becomes so ex- 
plosive that a dynamical event, such as the ejection of the 
helium envelope, may occur. The latter phenomenon would 
have interesting consequences for the binary systems of our 
interest. The ejection of the helium layer would leave a bare 
CO core as in the case of a spiral-in phase (Sect. HOt . yield- 
ing a lower pre-SN mass. 

Fig. |S] shows that both the maximum and central tem- 
perature of the 3.2 Mq model are already close to that re- 
quired for neon burning (panel (d)). The shell with maxi- 
mum temperature tends to move inward, and most proba- 
bly will reach the centre. The density at the neon burning 
front is lower than 10 8 gcm -3 (panel (6)). Hence, a 3.2 Mq 
helium star is expected to undergo core collapse due to pho- 
todisintegration, and most probably a dynamical event will 
not occur. The future of the 2.8 Mq model is less clear. 
The shell with maximum temperature moves outward which 
might leave part of the ONeMg core unburned. The central 
density is close to 10 8 gcm~ 3 . Whether a dynamical event 
occurs or not depends on how close the neon burning front 
is from the centre. Nomoto & Hashimoto (1988) argue that 
in a 2.8 Mq helium core, as a result of electron captures 
in the neon burning shell, the mean molecular weight fi c 
would increase above 2 which implies a lower Chandrasekhar 
mass. Therefore the gravothermal specific heat of the core 
remains negative and the central temperature will continue 
to increase. Neon burning will reach the centre and photo- 
disintegration will trigger the core collapse. 



3.3 Spiral-in phase vs. supernova explosion 

Helium stars more massive than 3.3 Mq do not undergo 
a CE and spiral-in phase in their late stage of evolution. 
The further evolution is straightforward. The core is weakly 
degenerate, and will collapse due to photodisintegration; the 
mass and period prior to the explosion are the same as those 
at the end of RLOF (i.e. M t = M Q , P t = P ). 

Helium stars of 2.8 - 3.3 Mq and those less massive 
than 3.8 Mq in very close orbits undergo a CE phase at 
the end of their evolution. The first question we will try to 
answer is whether the system will survive CE evolution, and 
how long the spiral-in phase will last. We ignore the possi- 
bility that the stars might experience a dynamical event due 
to the high density at the neon burning shell. Neon burn- 
ing (or the increase of temperature close to that required 
for neon ignition) and the high L v suggest that the helium 
star is already close to core collapse. The next question is 
whether the SN explosion takes place before or after the 
neutron star completes the spiraling-in process in the en- 
velope of the helium star. This depends on the competition 
between the spiral-in timescale and the remaining time until 
the explosion. 



3. 3. 1 The decay timescale of the spiral-in phase 

For a CE situation where the accretor is significantly less 
massive than the donor, so that the accretor does not cause 
a considerable perturbation to the structure of the donor, 
the orbital decay can be expressed as the change in orbital 
energy due to the drag force (Bondi & Hoyle 1944, Shima 
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et al. 1985), i.e. 
G M a Mns da 



(3) 



2 a 2 dt 

where M a is the mass in the giant interior to radius a, p a 
is the local density at separation a, and V is the relative 
velocity between the secondary and the CE. p is the Mach 
number, i.e. p — V/V B where V B is the speed of sound. The 
accretion radius Ra is approximated by 

2 G A/ns 



Ra = 



(4) 



V 2 + V 2 

The function £(/z) determines the dissipation rate, and is 
of order 2 - 4 in the supersonic case (Bondi & Hoyle 1944, 
Shima et al. 1985). We applied £(/i) = 2.5 which is the value 
for a steady-state set up after a violent perturbation (Bondi 
& Hoyle 1944). 

In order to estimate the timescale of the CE and spiral- 
in phase, we define the decay timescale as Tdccay = a/a and 
calculate the ratio of the decay timescale to the Keplerian 
timescale which is expressed as (Livio & Soker 1988 and 
references therein): 

Tdccay 



TKep 
1 

127T 



G( M ) 



M a + Mns 



Mns 



Vs 
Vkep 



(5) 



where Vkop is the Keplerian orbital velocity, and p a is the 
average density in the donor interior to radius a. G(p) is a 
function of the Mach number p given by 

1 (p 2 + lf 



G(p) 



(6) 



m M 3 

The parameter /3ce measures the importance of local (three- 
dimensional) effects in the spiraling-in process. If /?ce < 1, 
energy is deposited locally and neither spherical nor cylin- 
drical symmetry can be assumed. 

As a result of the deposition of orbital angular momen- 
tum into the envelope, the envelope is spun up, and the rel- 
ative velocity between the neutron star and the envelope is 
reduced. The drag force decreases, prolonging the spiraling- 
in process. We calculate the ratio of the spin-up timescale of 
the envelope to the decay timescale as (Livio & Soker 



7CE 



where 



Tdccay 



1.2 pf 



M a + Mns 



Mns 



Vb 
Vkep 



(7) 



Pa = -4 / r 4 p(r)dr (8) 

-Rin is the radius at the giant's core-envelope boundary. A 
considerable spin-up of the envelope is expected to occur if 
7CE < 1. 

The above estimation of the decay timescale is only 
valid if the donor is much more massive than the neutron 
star (cf. e.g. Iben & Livio 1993). In our calculation we have 
Mns ~ M , and hence the above derivation is not com- 
pletely valid, but nevertheless we use it to obtain a rough 
estimate of the decay timescale. 

Fig. E| presents the /3ce and 7ce parameters as a func- 
tion of radius in the envelope of a 2.8 Mq helium star with 
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Figure 6. The physical parameters in the common envelope of a 
2.8 Mq helium star with a 1.4 Mq neutron-star companion at the 
last model of our calculations for three different initial periods: 
the decay timescale with the scale on the right y-axis (dashed-), 
and the parameters /?ce (dotted-) and 7ce (solid-line). 



three different initial periods in the final models of our cal- 
culations. Here the relative velocity between the neutron 
star and the envelope of the helium star is assumed to be 
Keplerian. In all cases we find 7ce < 1, suggesting that 
spin-up of the envelope takes place at relatively early stages, 
and therefore the spiraling-in phase proceeds on a timescale 
longer than the estimated decay timescale. The 7ce param- 
eter increases and local effects become more important (/3ce 
decreases) with increasing initial period. This can be under- 
stood because a helium star in a wider system has more or 
less the same core mass as that in a closer system, but with 
a more extended envelope and higher density gradient, and 
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therefore is more centrally condensed (smaller p a but larger 

Pa/pa)- 

By integrating eq. ® over the envelope of the helium 
star, we estimate the timescale of the spiral-in phase, tdecay, 
as presented in Fig. [S] Note that idccay / Tdecay as defined 
above; Td eca y is the local decay timescale while idccay repre- 
sents the time required to spiral down to radius r. Taking 
these estimates at face value suggests that the spiraling-in 
process in the system with Pi = 10 d will last for 0.16 yr, 
while in the system with Pi = 2 d it lasts for 3 yr. However, 
because 7ce < 1 for all systems, in particular for P = 2 d , 
the actual spiral-in timescale will be longer. 

Furthermore, we have assumed that the structure of the 
envelope does not change during the spiral-in. If we consider 
that the deposition of orbital energy actually causes the en- 
velope to expand, then the expansion causes a decrease in 
Pa/ pa- This will lead to a slower orbital decay, by eq. ©. 
The quantity p a , which depends only on the structure of 
the envelope, decreases more significantly than p a (which 
is determined also by the core density) as the result of this 
expansion. Accordingly, 7ce decreases, and the spiral-in pro- 
ceeds even more slowly. Therefore, the decay timescale we 
have estimated above should be regarded as the lower limit. 
Note that Podsiadlowski (2001) has included the envelope 
expansion in his CE calculation and found that the spiral-in 
is initially very rapid and slows down after significant enve- 
lope expansion has taken place. 

3.3.2 The time left until the supernova explosion 

We first consider the possibility that the neon-burning shell 
does not reach the centre and therefore core collapse is trig- 
gered by electron captures. Because the core loses energy 
mainly in neutrinos, we can estimate the remaining time un- 
til the explosion by comparing the change in binding energy 
to the neutrino luminosity, L v . We calculate the binding en- 
ergy of the core of the last model, E , and estimate the bind- 
ing energy at the onset of electron captures, Ef. By assuming 
that at this point the core is completely degenerate, we solve 
Chandrasekhar's differential equation for the structure of a 
white dwarf (Kippenhahn & Weigert 1994). We applied p c 
= 3.7 x 10 9 g cm -3 as the central density at the onset of 
electron capture (Nomoto 1987, Miyaji et al. 1980, Miyaji 
& Nomoto 1987). By assuming that L„ remains constant af- 
ter the last model, we estimate the time until the explosion 
as AE/L V , where AE = E{ — E . Considering that L,, prob- 
ably increases instead of remaining constant, the timescale 
we derive gives an upper limit. For the 2.8 Mq helium star, 
we find that the maximum time is 70 yr for the system with 
P = 10 d and 103 yr for the system with P = 2 d . For a 
3.2 Mq helium star, the maximum time is 28 yr. 

Next we consider the possibility that the neon burning 
shell does reach the centre, and core collapse is triggered by 
photodisintegration which seems more likely than electron 
capture given the discussion in Sec. 13.21 A comparison with 
a detailed evolution calculation (Heger 2002, private com- 
munication) yields that a 1.696 Mq CO core will undergo 
core collapse ~ 20 yr after off-centre neon ignition. This core 
mass is about the same as the CO core mass of our 3.2 Mq 
model. Our 2.8 Mq helium star has a CO core of 1.45 Mq, 
and probably needs a longer time before undergoing core col- 
lapse. We conclude, therefore, that the 3.2 Mq helium star 



will undergo core collapse in 20 - 30 yr, and the 2.8 Mq 
model needs 20 - 100 yr before it collapses. 

Although we have tried to estimate the timescales for 
the spiral-in phase and for the star to undergo core collapse, 
the exact timescales remain uncertain. Therefore, there are 
two open possibilities for the outcome. If the orbital decay 
timescale is shorter than the time until collapse, the he- 
lium star explodes after the spiral-in phase is terminated. 
If the time until collapse is shorter than the orbital decay 
timescale, the SN explosion takes place before the neutron 
star completes the spiraling-in process (i.e. inside the CE). 
This has important consequences for the pre-SN mass and 
period of the system. 

3.3.3 The possible remants of the lower mass helium stars 

Assuming that there is enough time to spiral in before the 
explosion, we will now investigate whether the system sur- 
vives the CE and spiral-in phase. With the energy equation 
for CE evolution (Webbink 1984, de Kool 1990) 

<H Mcorc MnS 1 ,„x 

~ Mo M NS +2Menv/(??CE Ar L ) 

we calculated the post-CE separation, a t . Here M corc and 
M onv are the masses of the core and envelope, respectively. 
r-L — Rl/o is the dimensionless Roche radius of the helium 
star, and a Q is the pre-CE separation. ?7ce, the so-called 
efficiency parameter of CE, is taken to be 1. The parameter 
A describing the binding energy of the envelope to the core 
- approximated as in Dewi & Tauris (2000) - is calculated 
from the bottom of the convective helium envelope taking 
into account the gravitational binding energy only (if the 
internal energy is also taken into account, A and therefore 
a t would be larger). We find A is in the range 0.05 - 0.12 
(see Table HJ. We assume that the binary will survive the 
CE and spiral-in phase if the CO core does not immediately 
fill its Roche lobe in its new orbit. 

Assuming that the neutron star completes the spiral-in 
phase, the post-CE (i.e pre-SN) mass and period are pre- 
sented in Table |2| We find that all systems have final sep- 
arations larger than the radius of the core (a t >3P C ore), 
suggesting that they all survive the CE and spiral-in phase. 
Even if we consider that the core could expand by a factor 
of 2 after the envelope is peeled off, in most cases it still 
would not fill its Roche lobe. However, we should realize 
that the product 7]ce A in eq. © depends on the details of 
CE phase which are very uncertain and therefore the results 
in Table |5] give only approximate final separations. We find 
that the pre-SN orbit has a period of ~ d 01. In such a 
tight orbit, the neutron star moves with a very large orbital 
velocity (~ 10 3 kms -1 ). This velocity is higher than a plau- 
sible kick velocity of a few hundreds kms . Therefore the 
effect of a supernova kick is not as strong as in the case of 
systems in wider orbits. However, although an asymmetric 
explosion does not change the separation significantly, it can 
increase the eccentricity significantly, which in turn reduces 
the merger timescale (Belczynski, Bulik & Kalogera, 2002). 
If we assume that the explosion is symmetric and that the 
SN remnant has a mass of 1.4 Mq, then the mass which 
leaves the system is less than half of the initial total mass, 
and therefore all systems will remain bound (Blaauw 1961). 

After the explosion of the helium star, which is assumed 
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Table 2. The remnants of case BB and case BC evolutions from 2.8 - 3.2 Mq helium stars after surviving CE and spiral-in phase: the 
initial mass and period, the pre-CE mass and period, the mass and radius of the core, the binding energy parameter, the post-CE period 
and separation, the post-SN separation and eccentricity, and the merger timescale. 
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mass in a certain range, i.e. 2.8 - 3.3 Mq and 3.3 - 3.8 Mq in 
P O rb<0 d 25. With very short merger timescales, these DNSs 
would increase the detection rate of gravitational-wave ra- 
diation. Another implication is that the merger would take 
place relatively close to the host galaxy. Since a merger of 
a compact binary has long been thought to be one of the 
sources of 7-ray bursts, the existence of tight-orbit DNSs 
has important consequences for the understanding of 7-ray 
burst progenitors. We will try to estimate the birthrate of 
this new population in a forthcoming paper (Dewi, Pols & 
van den Heuvel, in preparation). 

The upper panel of Fig. |7| shows the masses and peri- 
ods at the end of our calculations (i.e. post-RLOF). Systems 
which go through a subsequent CE phase are indicated by 
open symbols. If the neutron star in those systems com- 
pletes the spiraling-in process, then their pre-SN (i.e. post- 
CE) masses and periods are presented in the lower panel. 
Since the explosion of the helium star can take place be- 
fore the neutron star completes the spiral-in phase, the up- 
per and lower panels of Fig. represent the maximum and 
minimum pre-SN mass and period, respectively, for systems 
which go through the CE phase. If the neutron star does 
not have enough time to complete the spiral-in before the 
helium star explodes, the pre-SN parameters are determined 



to leave a neutron-star remnant (with mass Mns,2), the or- 
bital evolution of the two neutron stars will be governed by 
gravitational-wave radiation as (Peters 1964) 

da _ 64 G 3 Mns M ns , 2 Mi 

dt ~ ~~5 ~ a 3(l _ e 2)7/2 



( l + R e * + E e *) ( io) 

V 24 96 J y ' 



de 304 G M NS M N s 2 M T (, 121 

e 



dt 15 c 5 a 4 (l-e 2 ) 5 / 2 

Assuming a symmetric explosion, we calculate the post- 
SN separation and eccentricity as will be explained in 
Sect. 14.21 These parameters are used to determine the 
merger timescales of these systems, which is calculated by 
integrating eqs. J10P and 111 U (for the complete equations 
please refer to Peters (1964)). This merger timescale is much 
shorter than that of the observed DNSs (see Table and 
also shorter than their characteristic age (i.e. ~ 300 Myr 
and 100 Myr, respectively, for B1913+16). This makes the 
probability for such systems to be observed very small. 

The existence of very tight-orbit DNSs has been also 
proposed by Belczynski et al. (2002), with the assumption 
that a CE phase occurs if the helium star is more massive 
than the neutron star. We have demonstrated in this work 
that a CE phase does occur, but only if the helium star has a 
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Figure 7. The masses and periods of the helium star-neutron 
star binaries at the end of our calculations. The upper panel gives 
the mass M and period P of the remnants before the spiral-in 
phase (after RLOF). The lower panel presents the mass Mt and 
period Pt after the spiral-in phase (prior to SN explosion). The 
star symbols represent the remnants of case BB mass transfer (see 
Paper I), and circles indicate the remnant of case BC evolution 
from this work. The solid symbols represent the remnants of he- 
lium stars which do not go through a CE phase, and the open 
ones indicate those which do. Note that for systems that do not 
go through a CE phase (solid symbols), Pt = Po and Mt = M . In 
the upper panel, thin lines connect the remnants of helium stars 
with the same initial mass Mj . The shaded area marks the region 
where a double neutron star can be produced by avoiding RLOF, 
taken from single helium stars calculation (Pols, in preparation) 
after taking into account the effect of stellar wind mass loss. 



by the situation prior to the core collapse, i.e. mass and 
period will be between the values given in the upper and 
lower panels. However, Podsiadlowski (2001) showed that 
the most rapid orbital decay takes place at the beginning of 
the spiral-in phase. Hence, although the core collapse occurs 
before the spiral-in is terminated, the pre-SN period is likely 
to be close to the minimum. Note that for systems that do 
not go through a CE phase (solid symbols), the mass and 
period in the upper panel are the same as in the lower panel. 

We have discussed in Sect. 13.21 the possibility of a dy- 
namical ejection of the helium layer due to explosive neon 
flashes for Mh c ~ 2.8 — 3.2 Mq. Although it is unclear 
whether this actually occurs, and we ignore it in the re- 
mainder of the paper, we will discuss briefly how such a dy- 
namical ejection might influence our results. If the ejection 



occurs before mass transfer becomes dynamically unstable, 
then we find the situation where the neutron star is orbiting 
the CO core, with mass as in the lower panel of Fig. |7] but 
with an orbit more like in the upper panel. If the ejection 
occurs after the spiral-in process is initiated, the future of 
the system depends on whether the ejection occurs after or 
before the completion of the spiral-in process. We will have 
the same situation as that where the core collapse occurs af- 
ter the completion of the spiral-in (i.e. the lower panel) if the 
dynamical ejection of the helium layer takes place after the 
neutron star terminates the spiral-in phase. If the ejection 
ensues before the completion of the spiral-in, then again the 
CO core remains (with a mass as in the lower panel) but the 
period will be between the upper and lower panels. 



3.3.4 The type of supernova explosion 

The final amount of helium left in the envelope probably 
determines whether the explosion will be observed as a type 
lb or a type Ic SN. The main observational criterion to 
distinguish between these types is the presence of helium 
in type lb and its absence in type Ic SN. The conclusion 
drawn in Paper I, that lower-mass helium stars and systems 
in close orbits are possible progenitors of type Ic SNs, and 
that higher-mass helium stars and systems in wide orbits 
produce type lb SN, is still valid here. 

We will discuss the possibility of the explosion type in 
the case of systems which undergo a CE and spiral-in phase. 
We consider the situation in which the SN explosion occurs 
after the neutron star completes the spiraling-in process, i.e. 
the whole helium envelope has been removed from the star. 
We discuss an extreme case, where the core collapses at the 
moment when the whole envelope is ejected. If the ejected 
matter is still surrounding the core, then the explosion will 
be observed as the helium-rich type lb SN. This situation 
resembles the case where the explosion occurs inside the CE. 

In the other extreme, we have the situation where the 
explosion occurs after the end of the CE phase. In the case 
of the 2.8 Mq model with Pi = 2 d , if we assume that the 
CE phase lasts for only a few years (Sect. l3~3.H and the ex- 
plosion takes place at most a hundred years after the start 
of the CE phase (Sect. |3~.3.2L then the core collapses less 
than a hundred years after the CE phase is terminated. We 
assume that the envelope is ejected with the escape velocity, 
i.e. 420 km s" 1 , and that during the SN explosion matter is 
ejected with a velocity of ~ 10 4 kms -1 . Using these veloc- 
ities, we find that in 5 years the SN shell will catch up the 
envelope matter, forming a shock front. The SN shell might 
interact with the envelope matter in a form of a ring around 
the SN remnant like in SN 1987A, which was observed a 
year after the discovery of the SN. The explosion itself may 
be very dim since the ejected mass is very small (less than 
~ 0.5 Mq if we assume that the SN remnant has a mass 
of 1.4 M Q , see Table|5J. The models of 2.8 Mq with larger 
orbital periods have slightly higher masses but much larger 
radii and therefore lower escape velocities, such that the 
timescales for the SN shell to catch up the envelope matter 
are shorter (~ 2 yr in Pi = 10 d ). In the models of 3.2 Mq, 
the escape velocities are large but the time left until the SN 
explosion is very short such that the interaction between the 
SN shell and the envelope matter takes place in ~ 1.5 yr. 
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Figure 8. Range of allowed mass and period of the immedi- 
ate pre-SN progenitor of double neutron-star binaries: B1913+16 
(dashed-), B1534+12 (dotted-), J1518+4904 (dash-dotted-), and 
J1811-1736 (solid-line). For each line style, upper and lower lines 
represent the maximum and minimum pre-SN periods, assuming 
that the ages of the pulsars are equal to their characteristic ages. 
The meaning of the symbols and shading are the same as in Fig. 171 
Diamonds indicate the pre-SN masses and periods for the case of 
symmetric supernovae for all four pulsars. 



4 THE FORMATION OF DOUBLE NEUTRON 
STAR BINARIES 

Following up on our attempt in Paper I, we will try to find 
constraints on the masses and separation of the progenitors 
of the observed galactic DNSs, and on the kick velocity that 
was imparted during the second SN explosion. This has been 
done before by various authors. However, as will be discussed 
later, here we reconsider their findings in the light of RLOF 
from helium stars which was not considered before. 

Table [3] lists the known galactic DNSs. Included in the 
table is the suggested DNS J1811-1736 (Lyne et al. 2000), 
although the possibility that the companion of the pulsar in 
this system is a main-sequence star or a red dwarf is still 
open (Mignani 2000). We plotted the allowed minimum and 
maximum pre-SN periods for each observed system in Fig.|H] 
which is expressed as (Flannery & van den Heuvel 1975) 



(1 - e f ) 3 M f P t 2 < M s P t 2 < (1 + eifMf P? 



(12) 



Here M s , Mf are the total pre- and post-SN masses; Pt, Pf 
are the pre- and post-SN periods. The values of Pf have been 
corrected for orbital decay during a time equal to the pulsar 
characteristic age (see Sect. l4~2l . Eq. Ill 2t comes from the 
assumption that the radius of the pre-SN orbit (the helium 
star-neutron star binary) must be between the periastron 
and apastron distance of the post-SN orbit. 



4.1 Previous studies 

Yamaoka, Shigeyama & Nomoto (1993) concluded that 
a symmetric explosion cannot explain the formation of 
B1534+12. Their conclusion was based on the argument that 
the pre-SN helium star mass is lower than 2.2 Mq, which 
is the critical mass for a helium star to become a neutron 
star (cf. Habets 1986). They argued that a kick velocity of 
300 - 460 kms" 1 (for B1913+16) and 160 - 260 kms" 1 (for 
B 1534+ 12) is required, assuming that the helium star does 



not fill its Roche lobe, which requires Mh b > 5 Mq in their 
adopted models . In the case where RLOF does occur they 
argue that the neutron star will spiral into the evelope of the 
helium star, such that the mass prior to SN explosion is the 
CO core mass which is smaller than the mass of the helium 
star. Since the required kick velocity increases with mass, 
a smaller kick velocity is needed for this case. Also in the 
situation where the helium star has experienced wind mass 
loss prior to the explosion, i.e. its final mass is much smaller 
than the initial helium star mass, a smaller kick velocity is 
required. 

Fryer & Kalogera (1997) also found that a symmetric 
explosion in the formation of the observed DNSs requires a 
pre-SN separation smaller than the maximum radius of the 
helium star, i.e. RLOF must have occured. They assumed 
this must result in a CE phase. They found that, if a neutron 
star is able to accrete matter above its Eddington limit in the 
CE phase, the timescale needed for a neutron star to collapse 
into black hole is much smaller than the timescale for a he- 
lium star to evolve from its maximum radius up to the explo- 
sion. They concluded that the neutron star will collapse into 
a black hole. Hence, in this case, a symmetric explosion fails 
to explain the existence of the DNSs. Minimum kick veloc- 
ities of 260, 220, and 50 kms" 1 (for B1913+16, B1534+12, 
and J1518+4904, respectively) are required, from a progen- 
itor with a t ~ 4.5 Rq, M t ~ 4.5 Mq (for B1913+16 and 
B1534+12) and a t ~ 30 Rq, M,~3M (for J1518+4904), 
in order to avoid RLOF from the helium star. As the com- 
pletion to this work, Wex, Kalogera & Kramer (2000) used 
the misalignment between the orbital angular momentum 
and the spin of B1913+16, together with its proper motion, 
to put more constraints on the kick velocity imparted during 
the second supernova. 

None of the above-mentioned works was based on de- 
tailed calculations of helium stars in binary systems. Fur- 
thermore, most of them were carried out assuming that the 
helium star does not fill its Roche lobe. Based on our cal- 
culations we will reinvestigate the kick velocity imparted 
at the birth of the pulsar's companion, allowing for RLOF 
from the helium stars. With this assumption, the low pre-SN 
mass found by e.g. Yamaoka et al. (1993) is not necessarily 
the initial mass of the helium star. Hence, we allow lower 
pre-SN masses. Also by allowing the helium star to fill its 
Roche lobe without causing the neutron star to collapse into 
a black hole, progenitors in closer orbits than those found 
by Fryer & Kalogera (1997) are possible. 

In Fig. |H] we compare the allowed range of pre-SN or- 
bital periods of the DNSs with the results of our model 
calculations, reproduced from Fig. |7| (lower panel). Fig. H 
shows that the possible progenitors of the short-orbit DNSs 
B1913+16 and B1534+12 are helium stars with M t > 
2.2 Mq, i.e. with initial mass Mh c > 3.3 Mq (solid sym- 
bols) , if indeed helium stars of lower mass leave DNSs with 
very small periods (open symbols). This means their progen- 
itors were main-sequence stars more massive than 12 Mq 
which underwent case B evolution, or stars more massive 
than 10 Mq which experienced case C mass transfer. 



4.2 Symmetric explosion 

In this section we reinvestigate whether a symmetric explo- 
sion can explain the observed DNSs, if the helium stars are 
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Table 3. The orbital parameters of the observed galactic double neutron-star pulsars: the masses of the pulsar and its companion, the 
total mass of the system in Mq, orbital period in days, eccentricity, and characteristic age in 10 8 yr. 



PSR 


M p 


M c 


Mf 


Pb 


e 




Reference 



1.08 Taylor & Weisberg 1989, 

Taylor 1992 
2.4 Wolszczan 1991, 

Stairs et al. 2002 
> 160 Nice, Sayer & Taylor 1996, 

Thorsett & Chakrabarty 1999 
9+2 Lyne et al. 2000 



B1913+16 1.4411 1.3874 2.82843 0.323 0.61713 

B1534+12 1.3332 1.3452 2.67843 0.421 0.27368 

J1518+4904 1.56t{5;^ L05 ±o:ii 2 - 62 8 - 634 0- 24 948 

J1811-1736 < 2.3 > 0.7 2.6 18.779 0.828 



allowed to undergo RLOF. By means of eqs. (1I0> and l|ll|l we 

calculate the evolution of eccentricity and separation of the 
DNSs B1913+16 and B1534+12 due to gravitational-wave 
radiation, back from the present time to twice their charac- 
teristic ages, as presented in Fig. as solid-line segments. 
The solid circle represents the eccentricity and separation 
for an age equal to the characteristic age. In a symmetric 
SN explosion, the post-SN eccentricity, ef, and separation, 
a,f, are related to the pre-SN mass, Mt, and separation, ot, 
as (Hills 1983) 

M t = et (M p + M c ) + M c (13) 
ot = 2 M c + M p - M t 

Of ~~ M p + M c [ ' 

where M p and M c are the masses of the pulsar and its com- 
panion, respectively. The star symbols in Fig. [5] represent 
the expected ef and Of for each of our model calculations, 
calculated from eqs. I13H and 1141 . 

The upper panel of Fig. [5] shows that in the case of a 
symmetric explosion, B1913+16 should have formed from a 
helium star with Mt = 3.13 — 3.42 Mq, depending on its 
true age (which is assumed to be less than twice its char- 
acteristic age). However, the separation is smaller than the 
minimum separation allowed for dynamically stable RLOF 
(cf. fig. 10 of Paper I). This implies that B1913+16 cannot be 
formed by assuming a symmetric explosion from helium star- 
neutron star binaries which undergo a mass-transfer phase, 
regardless of the exact age of the system. 

The eccentricity of B 1534+ 12 corresponds to M t = 
2.07 — 2.12 Mq. The eccentricity and separation of this pul- 
sar lie close to the remnants of our 3.4 Mq model (indi- 
cated by an arrow in the lower panel of Fig.|3J. However, we 
have discussed in Sect. 13. ll that the future of this particular 
mass is not clear. These systems may also undergo a CE and 
spiral-in phase and move to the lower-left part of the plane 
as in the case of lower-mass helium stars; and in that case 
cannot be considered as a possible progenitor of B 1534+1 2. 
Hence, a conclusion that a symmetric explosion can explain 
the formation of B1534+12 is marginal. Our results confirm 
previous work that it is more likely that both B1913+16 
and B1534+12 are formed by an asymmetric SN explosion, 
although the constraints on the kick velocity become much 
weaker if we allow for RLOF from the helium star - as will 
be discussed in Sect. 14.31 

The pre-SN mass and period (Mt, Pt) required in a 



symmetric explosion are indicated by diamond symbols in 
Fig. |S] for each of the four DNSs, assuming an age equal 
to the characteristic age. We find Mt = 1.96 Mq as the 
pre-SN mass of the wide DNS J1518+4904 if we assume a 
symmetric explosion. This mass (with a t = Of(l — ef)) is 
located in the region where helium stars of ~ 2.8 - 2.9 Mq 
will undergo a CE and spiral-in phase (compare Figs. and 
|SJ. Assuming that they complete the spiraling-in process, 
their position in the pre-SN mass-period plane will move 
to the lower-left part of Fig. |5] Even if they undergo a SN 
before completing the spiral-in, their final orbits are likely 
to be much closer than those at the onset of the CE phase. 
We conclude that a symmetric explosion cannot explain the 
formation of J1518+4904 and that it can only be produced 
from a helium star more massive than 2.8 Mq which avoids 
RLOF (shaded area in Fig.|8J, i.e. main-sequence stars more 
massive than 10 Mq in relatively wide orbits. 

Accurate masses of the components of J1811-1736 are 
not known yet. For our calculations we assume that the com- 
ponents have equal masses. With a symmetric explosion, we 
find Mt = 3.45 Mq, implying that the progenitor of J1811- 
1736 is a helium star of about 4.0 Mq. Helium stars with 
initial masses ~ 3.3 - 4.0 Mq in wide orbits which undergo 
marginal RLOF can be the progenitors of J1811-1736. How- 
ever, because of the much larger parameter space, a progen- 
itor that avoids RLOF is more likely. Whether a symmetric 
explosion can explain the formation of J1811-1736 depends 
on the exact masses of the components which in turn deter- 
mine the pre-SN at and M t . 



4.3 Asymmetric explosion 

If the kick velocity Vk makes an angle 9 with respect to 
the pre-explosion orbital velocity Vt, then we can write the 
relation between the pre-SN and post-SN separations, a t and 
a t , as (Hills 1983) 

£1 =2-^ ri + ^ 2 +2^cos0l (15) 
Of M f L J 

where M a = M t + M p , Mf = M c + M p , v = V k /V t and 
Vt = G M s /at. If we rotate Vk around Vt in a cone of apex 
angle 8 as is shown in Fig. 1101 then (f> is the location of Vk on 
this cone, such that <f> = corresponds to Vk in the original 
orbital plane and has a Cartesian component pointing radi- 
ally outward from the focus of the orbit (Hills 1983). The 
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Figure 9. The evolution of eccentricity and separation of 
B1913+16 and B1534+12 (solid lines) back from the present time 
to twice their characteristic ages (the values at the age equal to 
the characteristic age are indicated by the solid circle in the inset 
of each panel). Also plotted here, in solid stars, are the post- 
SN eccentricity and separation resulting from our calculations, 
assuming a symmetric supernova explosion. The pre-SN mass is 
indicated at the top x-axis. 
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Figure 10. The orientation of the kick velocity relative to the 
original orbital velocity Vt. 



larger than a critical angle, i.e. 

"2M f - M B (v 2 + 1) 



> Per = COS 



2vM s 



(17) 



Since cos 9 CI > — 1 an absolute maximum kick velocity can 
be derived from the above equation as 



VW = (1 + y , 2M i /M s ) V t 



(18) 



A kick with this magnitude has to be directed opposite to 
the orbital motion (8 = 180°). For M s > 2M f there is a 
minimum kick velocity 



= (1 - y/2 M t /M B ) V t 



(19) 



which also has to be directed at 6 — 180°. A 14, mm = 
requires M s < 2 Mi as derived by Blaauw (1961) for a sym- 
metric explosion. 

For a given observed DNS, we can find the pre-SN 
parameters by means of eqs. I15|l and (1 1 fit for a certain 
kick velocity (Vk,0,0). The kick direction is constrained by 
< sin 2 < 1. Independent of the magnitude of the kick 
velocity, the limit sin 2 = 1 gives eq. 1121 which can also be 
written as 



Of fi- 



ef. 



< a t < at (1 + ef) 



(20) 



These lower and upper limits in separation are presented as 
the thick horizontal lines in Fig. 1111 for the two closest DNSs. 
The limit sin 2 = for different kick velocity magnitudes is 
given by the thin lines. For a given kick velocity, the possible 
pre-SN parameters lie between the lines of sin 2 </> = 1 and 
that of sin 2 = (i.e. inside the shaded area marked by a 
certain kick magnitude). 



eccentricity of the post-SN orbit is given by 

Ga t Mt (1-e 2 ) =a 2 [V k 2 sin 2 (9 sin 2 + (V k cos (9 + V t f] 

or 

1 - e 2 = — — [1 + 2^cos6i + !/ 2 (cos 2 6i + sin 2 0sinV)l (16) 
at Mt 

After the explosion, a binary will remain bound if the right- 
hand side of eq. 11511 is positive. Hence, after an asymmetric 
explosion with a kick velocity Vk, a binary with pre-SN pa- 
rameters (M s , at) will not be disrupted if the angle 6 is 



4-3.1 Formation of DNS without a mass transfer phase 

We will first revisit the investigation of the formation of 
the DNSs assuming that the progenitors do not experience 
RLOF. A helium star will not fill its Roche lobe if the Roche 
radius is larger than its maximum radius, which defines a 
critical orbital separation a max - We plot this separation as 
the dashed line in Fig. 1111 This maximum separation is de- 
rived from the maximum radius taking into account wind 
mass loss as defined in eq. (|5J . This represents the situation 
where the helium star is formed after case B mass transfer, 
and has been losing mass by stellar winds during its evo- 
lution prior to RLOF. A minimum separation, below which 
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a helium star-neutron star binary cannot be formed, is de- 
fined by equating the Roche radius to the helium zero-age 
main sequence radius. Due to the wind mass loss, this min- 
imum separation increases according to eq. Q with a = 1 
and j3 = 0. We plot the region with separation less than the 
minimum separation at the end of the evolution of a helium 
star against its final mass as the hatched area. Helium star- 
neutron star binaries cannot be formed in this part of the 
diagram. These constraints leave the region above the maxi- 
mum (dashed line) and minimum (hatched area) separations 
as the allowed pre-SN parameters in Fig. II II (however, these 
constraints depend very much on the choice of the wind 
mass-loss rate). 

Without RLOF from the helium stars, we find the same 
minimum kick velocity as Fryer & Kalogera (1997), i.e. 260 
and 220 kms -1 for B1913+16 and B1534+12, respectively 
(see the shaded area above the dashed line in Fig. 1111 . The 
kick must be directed backwards, with 9 > 130°. For a low 
14 (< 425 kms" 1 for B1913+16 and < 275 kms -1 for 
B 1534+ 12) a kick in any azimuthal direction <j> is allowed 
to form the binaries with the observed parameters, but the 
higher Vk the more restricted the allowed orientation of the 
kick. 

The inclusion of the minimum separation (hatched area 
in Fig. II IB as an additional constraint, which has not been 
taken into account in previous works, enables us to deter- 
mine the maximum pre-SN mass. The maximum pre-SN sep- 
aration of B1913+16 intersects the line of ot = a m in at M% ~ 
5.8 Mq, which is the final mass of a ~ 16 Mq helium star 
after stellar wind mass loss. This means that B1913+16 can- 
not have formed from a helium star initially more massive 
than 16 Mq. Similarly, we derive ~ 13 Mq as the upper 
limit for the helium star progenitor of B1534+12. The max- 
imum pre-SN mass also implies a maximum kick velocity, 
from eq. (I18L which is 1 230 kms" 1 for both systems. 

In discussing the above situation, we consider that the 
helium star-neutron star binaries are produced at the time 
when the helium star is on its zero-age main sequence, i.e. 
the remnant of a Be/X-ray binary which went through a CE 
phase initiated in case B. If the CE phase is initiated in case 
C (i.e. after the termination of helium core burning), the he- 
lium core does not experience significant wind mass loss and 
therefore has a maximum radius equal to that of a helium 
star evolving without wind mass loss. This maximum radius 
is plotted as the dashed line in Fig. 1121 For case C remnants 
there is no minimum separation which is constrained by the 
helium zero-age main sequence radius. Whether RLOF can 
occur or not depends on the situation at the cessation of the 
CE phase which is very uncertain. We will consider the case 
where the case C remnant does not fill its Roche lobe (the 
area above the dashed line in Fig. 1121 . This area is located 
outside the range of possible pre-SN periods of B1534+12. 
Therefore we conclude that B1534+12 cannot be formed by 
case C mass transfer from Be/X-ray binaries without a fur- 
ther RLOF phase. To produce B1913+16 from this type of 
remnant, a very high 9 > 150° is needed. The upper limit 
to the pre-SN mass, and the maximum kick velocity, are 
restricted by the threshold mass for black hole formation, 
which is still uncertain. 

For the two wide-orbit DNSs, J1518+4904 and J1811- 
1736, the maximum separation for RLOF for case C rem- 
nants coincides with case B, because helium stars of such 



mass (2.5<M 4 /Mq <4.0, see Fig. lose very little mass 
in a wind. Hence, the minimum kick velocities are the same 
in both cases, i.e. 50 kms" 1 for J1518+4904, as also found 
by Fryer & Kalogera (1997), and 10 kms" 1 for J1811-1736 
(the latter depending on the uncertain component masses for 
this system). The maximum kick velocity, again, depends on 
the threshold mass for black hole formation. 

4.3.2 Formation of DNS with RLOF 

The region in the (Mt, a t ) plane where close-orbit DNSs can 
be formed through RLOF is limited, as shown in Fig. 1111 It 
is bounded at small separations because mass transfer be- 
comes dynamically unstable leading to a merger (below the 
dotted line). In Paper I we found that mass transfer from 
Mec > 6.7 Mq is dynamically unstable leading to a CE 
phase. This range of mass is shown as the region to the right 
of the dotted line (with M t >4 Mq). We have not carried 
out calculations for this range of mass. However, we suggest 
that such a CE phase would lead to a merger or a formation 
of DNS in tight orbit, depending on the binding energy of 
the envelope. At low Mt, to the left of the dash-dotted line, 
systems experience a CE phase after RLOF. This leaves the 
region between the dotted, dashed, and dash-dotted line, 
in which the results of our RLOF calculations lie (the star 
symbols). We do not present the star symbols in Fig. 1121 
because the calculations we have done are the remnants of 
case B mass transfer. Although we have not computed the 
evolution of the remnants of case C mass transfer, we ex- 
pect that the limits for dynamically unstable mass transfer 
(dotted line) and for systems undergoing a CE phase after 
RLOF (dashed-dotted line) are the same as for the remnants 
of case B evolution. 

We find that B1913+16 can be formed by an asym- 
metric explosion with minimum kick velocity of 70 kms" 1 , 
which must have been directed along the orbital motion 
(9 < 20°). A kick velocity as low as 10 kms -1 is enough 
to produce B1534+12, which requires 9 > 85°. A large Vk 
requires 9 close to 180°, and can only be imparted in a very 
close orbit (at ~ Of (1 — 6f), cf. eq. 1151 1. With this a t and 
using eq. 1181 we find 1810 kms" 1 and 1310 kms" 1 as the 
maximum kick velocities that can produce B1913+16 and 
B1534+12, respectively. To produce a system with post-SN 
parameters (Mf, at) from a binary with pre-SN parameters 
(M s , ot), it can be seen from eq. I|160 that with the same 
kick velocity, a decreases with 9. Since B1534+12 has al- 
most the same mass and separation as B1913+16 but lower 
eccentricity, in general, with the same kick velocity, a larger 
angle 9 is needed to produce B 1534+ 12 than that to form 
B1913+16. 

We have discussed in Sect. l3.3"21 that there are two open 
possibilities for the fate of helium stars which go through a 
CE phase. If the core collapses after the spiral-in phase is 
completed, then the system becomes a double neutron star 
in a tight orbit, i.e. the situation is as described in the previ- 
ous paragraph. If the helium star explodes while the neutron 
star is still spiraling-in, the pre-SN mass and separation will 
be somewhere between M and Mt and between P and P t 
(cf. the upper and lower panels of Fig. [7J, i.e. the region 
to the left of the dash-dotted line in Fig. 1111 may not be 
empty. The inclusion of this area with lower mass would 
allow a wider range of kick velocities and directions. The 
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Figure 11. The possible magnitude of the kick velocity imparted 
at the birth of the second neutron star in B1913+16 (upper) 
and B1534+12 (lower panel), with the assumption that the pul- 
sar age is the same as its characteristic age. The star symbols 
represent the pre-SN parameters of our case BB and BC calcu- 
lations which do not go through a CE phase. Thick horizontal 
lines mark the minimum and maximum pre-SN separations con- 
strained by sin 2 = 1. Thin lines are obtained from the constraint 
of sin 2 = 0. The solid circle (indicated by an arrow) gives the 
pre-SN parameter if the explosion is symmetric (V^ = 0). Below 
and to the right of the dotted line, mass transfer occurs on the 
dynamical timescale leading to a merger. To the left of the dash- 
dotted line is the region where mass transfer ends in a CE phase. 
The dashed line represents the maximum separation derived from 
the maximum radius of a helium star produced from case B mass 
transfer, after taking into account the stellar wind mass loss. The 
hatched area indicates the region where the separation is smaller 
than the minimum separation derived from the helium zero-age 
main sequence. The corresponding kick velocities for each panels 
are presented at the bottom of the figure. 
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Figure 12. Same as Fig. 1111 The dashed line now represents the 
maximum separation derived from the maximum radius of helium 
star produced from case C mass transfer, i.e. without stellar wind 
mass loss. 

magnitude and orientation of the minimum kick velocity for 
B1913+16 remains the same. In this situation B1534+12 
can be produced by a symmetric explosion or with very low 
kick velocity in all directions 6 and (j>. 

Comparing our results with those of Yamaoka et al. 
(1993), Fryer & Kalogera (1997), Wex et al. (2000), and 
Sect. 14.3.11 we can see that by allowing stable RLOF from 
the helium star, systems with lower pre-SN mass and lower 
pre-SN separation than derived in previous works can also be 
possible progenitors of B1913+16 and B1534+12. We have 
shown that a low kick velocity can then explain the forma- 
tion of these DNSs. 



5 CONCLUSIONS 

We have evolved helium stars with masses in the interval 2.8 
- 6.4 Mq with a 1.4 Mq neutron-star companion, in which 
the helium star fills its Roche lobe during carbon core burn- 
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ing or beyond (case BC mass transfer). This is the comple- 
tion to our earlier work, i.e. mass transfer during helium core 
burning (case BA) and helium shell burning (case BB evo- 
lution) . We studied the late stage of evolution of the helium 
star-neutron star binaries as well as the possible remnants 
of the systems. 

Case BB and BC mass transfer in helium stars of 2.8 
- 3.3 Mq as well as from 3.3 < Mhc/Mq < 3.8 in very 
close orbits (P<0 d 25), will end up in a CE phase towards 
the end of their evolution, just before the expected SN ex- 
plosion. These systems originate from main-sequence stars 
with masses of 10 - 12 Mq which underwent case B evolu- 
tion, or 9 - 10 Mq which experienced case C mass transfer. 
We found that all systems are able to survive the spiral-in 
phase, producing very tight DNSs with P ~ d 01. These sys- 
tems probably will not be observed due to their very short 
merger timescale, and would have important consequences 
for the detection rate of gravitational-wave radiation and for 
the understanding of 7-ray burst progenitors. On the other 
hand, there is a possibility that the helium star will explode 
before the neutron star completes the spiraling-in process, 
resulting in a SN explosion inside a CE. 

For A/ Ho > 3.3 M or M Ho > 3.8 M in close or- 
bits, we conclude that a CE and spiral-in phase does not 
occur. These systems will produce DNSs with periods of 
Orl — 1 , which suggests they are candidate progenitors of 
B1913+16 and B1534+12. The pre-SN mass is larger than 
2.2 Mq. These systems originate from main-sequence stars 
more massive than 12 Mq which underwent case B evolu- 
tion, or more massive than 10 Mq which experienced case 
C mass transfer. 

We have also studied the second SN explosion to in- 
vestigate whether a kick velocity is required at the birth of 
the young neutron star. DNS B1913+16 cannot be formed 
by a symmetric explosion. A minimum kick of 70 kms -1 
in the direction of the orbital velocity is needed to produce 
the observed orbital parameters. A symmetric explosion can 
explain the formation of B1534+12 only if the progenitor, a 
3.4 Mq helium star, avoids a CE and spiral-in phase. Our 
calculations are inconclusive on this issue. A higher mass 
progenitor needs a minimum kick velocity of 10 kms -1 in 
the direction opposite to the orbital velocity. 

The wide-orbit DNS J1518+4904 (and probably J1811- 
1736) can only be produced from helium star-neutron star 
systems which did not go through a mass-transfer phase; 
with the helium stars more massive than 2.5 Mq, i.e. main- 
sequence stars more massive than 10 Mq in relatively wide 
orbit. Without a mass-transfer phase, an asymmetric explo- 
sion with a minimum kick velocity of 50 and 10 kms -1 (for 
J1518+4904 and J1811-1736) is needed. 
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